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Any non-gravitational coupling between dark matter and dark energy modifies the cosmological dynamics. 
Interactions in the dark sector are considered to be relevant to address the coincidence problem. Moreover, 
in various models the observed accelerated expansion of the Universe is a pure interaction phenomenon. Here 
we review recent approaches in which a coupling between both dark components is crucial for the evolution 
of the Universe. 
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I. INTRODUCTION 

According to the prevailing view, our present Universe 
is dynamically dominated two so far unknown compo- 
nents, dark matter (DM) with an energy density p m and 
dark energy (DE) with a density p x . DM is usually mod- 
eled as a pressureless fluid, while DE is characterized 
by an equation of state (EoS) parameter w. The de- 
termination of this parameters is one of the major chal- 
lenges in the field. Current observations are consistent 
with a spatially flat universe with fractions of 72% DE 
and 28% matter (including DM and baryons). The pre- 
ferred model is the ACDM model (A denotes the cos- 
mological constant, equivalent to a constant energy den- 
sity p x = pa and "CDM" stands for cold dark mat- 
ter) for which w = — 1. Because of the cosmological 
constant problem in its different facets, competing ap- 
proaches were developed which either "dynamize" the 
cosmological constant or modify General Relativity. Here 
we review various approaches in which DM and DE do 
not evolve separately but interact with each other non- 
gravitationally. This gives rise to a richer cosmological 
dynamics compared with non-interacting models. More- 
over, disregarding the potential existence of a coupling 
between both dark components may result in a misled 
interpretation of observational data 1 . 



II. THE DARK SECTOR OF THE UNIVERSE 
A. The ACDM model 

The basic equations of the ACDM model are 
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and 
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where if = - is the Hubble rate and a is the scale factor 

a 

of the Robertson- Walker metric. 



Introducing the critical density p c = ||^? and the frac- 
tional quantities 
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the Friedmann equation (JT|) can be rewritten as 
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The current data are consistent with (the subindex de- 
notes present values) 



fi ra o ~ 0.28, f^Ao ~ 0.72 
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With the definitions A = 8ir G pa and pa = — Pa, the cos- 
mological constant can be seen as a "fluid" with negative 
pressure. Then 



8ttG 



k 



4ttG 
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where p — p m + pa and p = p m + p\. On the other 
hand, the cosmological constant is dynamically equiv- 
alent to the lowest energy state of a scalar field with 
energy-momentum tensor 
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Namely, considering the configuration 

it becomes obvious, that in the field equations 



R lk ~^g lk R + Ag„ 
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there appears an effective cosmological "constant" 
Aeff — A + 8irG p V ac- The density p vac of the vacuum 
energy (recall E vac — ^ Hu) and u) = \Jk 2 + to 2 ) is 
calculated via 
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where we have assumed % = 1 and m — 0. For a cutoff 
fe c at the Planck scale one finds p vac ~ 10 76 (GeV) 4 s» 
10 114 ^^-. The observed value for the effective cosmologi- 
cal "constant" is ,t ~ 10 _47,/ ' 1 



10" 



9iES . This 

cm 3 



is p Aeff — iir''(GeV) 
large discrepancy of 123 orders of magnitude constitutes 
the cosmological constant problem. 

Another problem results from the circumstance that 
the ratio of the energy densities, r = £ffiL , has a present 

value of ro ~ | which is of the order of unity. If one 
takes into account that, for the ACDM model, p x = p\ = 
const and p m oc a~ 3 , it follows that r oc a~ 3 . A value of 
ro of the order of one seems to single out a very special 
period. In general, r differs from unity by many orders 
of magnitude. This is called the coincidence problem. 



III. INTERACTING DARK ENERGY 

In the ACDM model and in most other approaches 
DM and DE are considered to be uncoupled components 
with separately conserved energy-momentum balances. 
However, Einstein's equations only require total energy- 
momentum conservation. This freedom can be used to 
construct alternative models of the dark sector in which 
DM and DE also interact non-gravitationally. We start 
with some general relations before considering specific 
models. We assume the total energy-momentum tensor 
to have a perfect fluid structure, = puiu^ +phik with 
T% = 0. Further, we assume a split into two compo- 
nents, 



we assume them to coincide: — u x — u a . Restricting 
ourselves for the moment to the background, the individ- 
ual energy balances are 

p m + 3Hp m = Q, p x + 3H(l+w)p x = -Q , (17) 

where Q = —u a Q a is considered to be a phenomenolog- 
ical quantity. Obviously, Q does not directly enter the 
Hubble rate and the deceleration parameter. Explicitly, 
it only appears in the second derivation of the Hubble 
rate^, 
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Its influence on the dynamics may be quantified by the 
"statefinder parameter" ("jerk") 
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This parameter enters the luminosity distance 

dz 



d L = (l + z) 



in third order in the redshift: 
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with perfect-fluid structures as well, i.e., Tjf = 



PAU" A u n A +PAh l A and h A = g' n + u" A u A , wtiere A — m, x. 
A coupling of both components is described by 



While the ACDM model has j = jo = 1, any interacting 
model has jo ^ 1 in general, which allows a discrimina- 
tion between models that share the same values of Ho e 
<7o- 
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A. Scaling cosmology 



which guarantees the overall energy-momentum conser- 
vation. The separate energy balance equations are 



and 



Pm,.aU a m + 6 m (p m + Pm) = ~U ma Q a (13) 



Px,aU a x + ®x {fix + Px) = U xa Q a , (14) 



where Qa = u A . k . Likewise, we have the momentum 
balances 



and 



(Pm "I" Pm) U m + Pm,ih m 



(px+Px)u a x +PxAx =-K i Q i 



(15) 
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Starting point is an ansatz for the ratio of the energy 
densities of DM and DE, 
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(22) 



with a phenomenological parameter £. The present value 
of the scale factor has been set to a = 1 . The ACDM 
model has £ = 3. A stationary ratio is characterized 
by £, = 0. The severity of the coincidence problem is 
quantified by £. According to Ref. 0, any £ < 3 alleviates 
the coincidence problem. It has been showr*^, that an 
interaction 



Q 
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l + r (l + z)^ 
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Generally, the four velocities of the components are dif- 
ferent from each other and from the total four velocity. 
For the homogeneous and isotropic background, however, 



can produce any scaling behavior of the type (|22[) . The 
interesting special case of a stationary ratio r = r = 
constant results in a power-law solution p Xl p m oc a~ u , 
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where v = 3 1 + r °+ w . This solution is able to describe 

l+r 

an accelerated expansion under the condition 3w < 
— (1 + 7'o). As any power-law solution, the dynamics can 
be described in terms of a scalar field with an exponen- 
tial potential 4 . The interesting feature here is that a spe- 
cific interaction can give rise to a stationary ratio of the 
energy densities which is consistent with an accelerated 
expansion of the Universe. 

For a toy model w = — 1, £ = 1 the impact of the 
interaction on observational quantities can be visualized 
as follows. The total energy density of this model is 



Pm + Px = P = - — — [l + r Q a x ] 
(1+ro) 



while one has 

P(KCDM) 



(1 + ro) 



[l + r a- 3 ] 



(24) 



(25) 



for the ACDM model. The sum of powers in (|23|) is re- 
placed here by the power of a sum in (pM)) . This difference 
affects the luminosity distance which, for our toy model 
becomes 
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in contrast to its ACDM counterpart 
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ro)z 2 ■ 
. (27) 

Both expressions differ in third order in the redshift. 

What one would like to have is a dynamical evolution 
toward a constant ratio r. It has been shown in Ref. [Bj 
that indeed a suitable interaction between DM and DE 
can drive the transition from a decelerated, matter dom- 
inated expansion to a DE dominated period with accel- 
erated expansion and a stable, stationary ratio of the 
energy densities. 



B. Holographic dark energy 

An interaction modifies the a~ 3 behavior of the matter 
energy density to 

p m = Pmoa~ 3 f (a) , (28) 

where the function f(a) encodes the influence of the in- 
teraction. The corresponding DE balance can be written 
in terms of an effective EoS parameter w e ^ , 

Px = -3H(l + w eff ) Px , w eff = w +-L-r . (29) 

In this context one can consider a constant energy density 
ratio: 



r = const w' 



eff 



f 



3Hf 



W 



eff 



1 + r 



(30) 



The effective EoS coincides with the total EoS - = w e f f . 

p 

Via Friedmann's equation 

3H 2 = 8ttG (p m + Px )=8nG (l + r)p x (31) 

it follows that for r = const we have necessarily p x oc H 2 . 
Now the question arises, whether there exist DE models 
with p x oc H 2 . This leads us to holographic DE models. 
Holographic DE models are characterized by an infrared 
cutoff length L which is related to an ultraviolet cutoff, 
corresponding to the DE density. As a consequence one 
obtains a dependence p x oc L~ 2 . The underlying idea 
here is that the energy in a volume L 3 should not exceed 
the energy of a black hole of the same size^: 



p x L 3 < 



L 



IvrG 



=> Px = 



3 c 2 



8ttG L 2 



(32) 



where c 2 is a constant of the order of unity. A choice 
L = H -1 then provides us with 



Px 



3 c 2 H 2 
8ttG 



(33) 



This expression was shown£ to yield a correct order of 
magnitude of the presently observed value of the cosmo- 
logical term. From (|3T)]) it follows that under this condi- 
tion w e ff and w are only different from zero for a non- 
vanishing interaction. In the non-interacting limit the 
only possible EoS is w — 0, which is incompatible with 
an accelerating universe. A negative EoS parameter is a 
pure interaction effect. This becomes particularly clear 
if we look at the expression 



(34) 



for the deceleration parameter q. Obviously, the sign of 




q depends on the ratio -gj. For y < H we have q > 0, 
while q < holds for y > H. The interesting ques- 
tion is, whether an evolution from y < H to ^ > H is 
possible. This would be equivalent to a transition from 
decelerated to accelerated expansion as a pure interac- 
tion phenomenon. Let us consider to this purpose the 
simple example of a constant interaction rate T where 

r 



-2-. In this case, the Hubble rate isi 



H = H 



3iZ r 



1 - 



3H r 



-3/2 



(35) 



which represents a special generalized Chaplygin gas. 
The evolution of the EoS parameter is shown in Fig. 
[TJ which confirms that there is indeed a transition from 
in » at high redshifts to an EoS parameter w ~ —1 at 
present. 

Introducing the density contrast S m — where p m 
is the first-order DM density perturbation and assuming 
for simplicity that the corresponding DE perturbations 
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S x = are related to the matter perturbations by S x = 
aS m with a constant a, we obtain the following closed 
equation for the DM density contrast 7 : 



x" i ^ 



r 2(1 -a) r 



3Hr 

a 



2a 2 



1 



1 - 



3 Hr 
4(1 -a) T 



Hr r 



S m = 



(36) 



The prime denotes derivation with respect to the scale 
factor. Alternatively, the fluctuation dynamics is charac- 
terized by the growth rate parameter / = ^fefr , which 
obeys the equation 



f' + f+l(l + ±-)j.~l = ». (87. 

In the limit T = we have / = 1, i.e., we recover the 
Einstein-de Sitter behavior S m oc a. For the model based 
on (|35|) and (|37p . the quantity / is visualized and com- 
pared with observational data in Fig. [2] 




FIG. 1. Evolution of the equation of state parameter of dark 
energy for the best fit model^. The red swath indicates the 
region obtained by including the la uncertainties of the con- 
strained parameters Q x and F/Hq- 



la ACDM 

Q x -0.745 r/H =0.563 

ACDM - 



FIG. 2. Growth function vs. redshift for the best fit holo- 
graphic model (solid line)£. Also shown is the prediction of 
the ACDM model (dashed line). 



C. Transient acceleration 

There are results in the literature which suggest that 
cosmic acceleration could be a transient phenomenon^. A 
dynamics of this type may be the result of an interaction 
between DM and DE as can be shown by the toy model 
considered in Ref. ITol . Let us parametrize the interaction 
according to 

f(a) = l+g (a) (38) 

and consider the special case 

w = -1 , g(a) = 7a 5 exp(-a 2 /ff 2 ) , (39) 

where 7 is an interaction constant. The DE density be- 
comes 



Px = P e jj - 7 ^ exp (-a 2 /- 2 ) ( 
1 + .90 \ 

with an effective cosmological constant 



a 2 -Vj (4<„ 



eff 
Hx 



Px Q - 2 7 r 



9o 



exp(-l/ < 7 2 ) 



(41) 



The interaction renormalizes the bare (interaction-free) 
value p Xo . A transient acceleration is only possible for 
p e JJ — since otherwise the constant would always pre- 
vail in the long-time limit. This means, for accelerated 
expansion to be a transient phenomenon, part of the in- 
teraction has to cancel the bare cosmological constant. 
Under this condition and with 



K 



SttG p mo 



3tf 2 



1 



9o 



we obtain 



1- §ifa 2 exp(-l/a 2 ) 



-Hi {-3 J FCexp(-a 2 /cr 2 ) [a 2 - a 2 ] } 



(42) 

(43) 
(44) 



This may be compared with the corresponding expression 
for the ACDM model 



1 - fi/ 



- 20a 



(45) 



Since one expects an alternative model not to deviate 
too strongly from the ACDM model at the present time, 
the comparison between (j4"4")l and (|4"5"j) suggests K to be 
positive and the interaction term to play a similar role 
as the cosmological constant A. In other words, the role 
of the interaction is twofold. As already mentioned, it 
has to cancel the bare cosmological constant. But at the 
same time it has to induce an accelerated expansion by 
itself. The expression (l4"4l guarantees an early (a <C 1) 
decelerated expansion for f2 mo > K exp(— 1/ct 2 ) which 
represents an upper limit on the interaction strength K. 
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FIG. 3. Acceleration parameter for the best fit scenario in units 
of H% (taken from Ref.[l(|. 



To reproduce the presently observed accelerated expan- 
sion requires the condition 



aH 2 



>0 ^ K exp(-l/cr 2 ) 



> 5 , (46) 



which represents a lower limit on the interaction strength. 
It follows that there is an admissible range 



2 e 1 /- 2 



< K < 



2e^ 2 

3<7 2 



(47) 



for K. From the data of the GOLD sample we found 
best-fit values a = 5.000, K = 0.019, h = 0.644, Q x0 = 
0.678 and Xmin — 1-129 which have to be contrasted with 
the corresponding values for the ACDM model which are 
h = 0.64, n x0 = n M = 0.65 and X 2 min = 1-128. Fig. [3] 
shows that this model indeed describes an early transition 
from decelerated to accelerated expansion together with 
a future transition back to decelerated expansion 10 . 



D. Decaying vacuum energy 

A still alternative scenario is based on a prescribed be- 
havior of the decay of the cosmological term, interpreted 
as vacuum energy. The only preferred time scale in a 
homogeneous and isotropic universe is the Hubble time 
H~ l . It is tempting to associate a supposed vacuum de- 
cay with this scale. The simplest case therefore is a de- 
pendence px oc H . This relation has also some support 
from QCI3ii, but our approach here is phenomenological. 
The dark-energy density is characterized by 



px 



-e 



Px 



-px 



e = 



(48) 



which, in the homogeneous and isotropic background re- 
duces to p x — aH, p x = —aH and = 3H. For the 
Hubble rate we find 



H — Hq 1 — £l m Q + ^ ra o« 



-3/2 



(49) 



where we used that a — jf- (1 — f2 m o). The ratio of the 
energy densities scales as 



Pm ^mO 

Px 1 — ^mO 



-3/2 



(50) 



To consider the perturbation dynamics we split the 
interaction term Q l in (fT2")) - (fT6"| into parts parallel and 
perpendicular to u\ 

Q l = u l Q + Q l , (51) 

such that 

Q = -u,Q l , Q l = hlQ a , Ul Q l = . (52) 
The basic set of equations then is 

Pro + ©Pro = - g © , PmUa = ~^K@,i (53) 

and the Raychaudhuri equation 

6 + \e 2 - u" a + 4ttG (p + 3p) = . (54) 

The line element for scalar perturbations is written as 

ds 2 = - (1 + 20) At 2 + 2a 2 F a dtdx a 

W [(1 - 2ip) 5 af) + 2E, aP ] dx a dx . (55) 

For the spatial components of the four- velocity we write 

a 2 u^ + a 2 F,„ = u^ = , (56) 

thus introducing the velocity potential v. The first-order 
matter energy balance becomes 

5 m + 4xGa5 m - (-3H + 4nGa) + 6 = — . (57) 

Pm 

Now it is convenient to formulate the perturbation dy- 
namics in terms of the gauge-invariant quantities 

9 C = e + Ov , 5 c m = S m + , 8% = Sx + P—v . (58) 

Pm Px 

Then the energy- and momentum balances are combined 
into 



fa + ^*ro + J& = , 



where 



aH ( AttG a aH k 2 
J= 1+ - 1- 



3p r , 



3 H 3p m a 2 H 2 



(59) 



(60) 



is a scale-dependent (k denotes the comoving wavenum- 
ber) quantity. Differentiation of (|59[) and combination 
with the first-order Raychaudhuri equation for C yields 
the second-order equation^ 2 - 



6% + g(a,k)6 c m + f(a,k)6 c m = 



(61) 



Interacting Dark Energy 



6 



s 

g 0.06 



0.0000 0.0005 0.OO10 0.0015 0.0020 

k (h/Mpc) 



FIG. 4. Relative power spectrum as a function of k on large scales 
for f2 m o = 0.3 (taken from Ref. 1 12ft . 



where the coefficients 17(0, k) and /(a, k) are known, 
scale-dependent background functions (see Ref. It 
is an advantage of this model that it allows us to calcu- 
late the dark-energy perturbations in terms of the matter 
perturbations and their first derivative. The relevant re- 
lation is 



8° = -+- 
3J 



a6Z 



1 - n 



M0 



1 - fijwo + &Moa 3/2 



(62) 



From the definition (|60p of J it follows that on small 

scales =j 3> H 2 one has \ J\ >• 1, equivalent to 5%. <C 5^, 
i.e., the DE perturbations are negligible on these scales. 

Pressure perturbations in interacting two-component 
models are generally non-adiabatic. The relevant quan- 
tity is 



The relative power-spectrum, based on equations (j6"Tj) 
and (|6"2"|) . is shown in Fig. |U It indicates an increase of 
DE perturbations on very large scales. 

IV. SUMMARY 



We have summarized here various aspects of the poten- 
tial role of interactions in the cosmological dark sector. 
In the simplest case this comprises third-order redshift 
corrections in the luminosity distances of supernovae of 
type la. But a suitable coupling may also replace the 
cosmological constant and generate a phase of acceler- 
ated expansion by itself. This can happen, e.g., in the 
context of holographic DE models and in models of tran- 
sient acceleration. In the latter case the future evolution 
of the Universe is qualitatively different from that pre- 
dicted by the ACDM model. A decaying-vacuum-energy 
approach is another scenario in which interactions have a 
noticeable impact on the cosmological dynamics. Quite 
generally, interacting models are characterized by non- 
adiabatic pressure perturbations, even though the indi- 
vidual components are adiabatic on their own. 
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For adiabatic perturbations p — jp = is valid 
At high redshift, equation (|5T|) reduces to 



2a 



2a 2 







(a « 1) , (64) 



which coincides with the corresponding equation for an 
Einstcin-de Sitter universe. The non-adiabaticity is very 
small for a <C 1 but it vanishes exactly only for f2 m Q = 1 . 
Under this condition it is possible to impose adiabatic 
initial conditions for the evolution of the matter pertur- 
bations which are then used to calculate the matter power 
spectrum. This corresponds to the circumstance that DE 
perturbation are very small for small values of the scale 
factor, 



\Px 



1 1 

3~ 



^m0 



,3/2 



^m0 



(«« 1) 
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